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Abstract 

Inclusive J /if} productions in e + e~ annihilation is studied in the framework 
of NRQCD. We first review the leading-order calculations of the cross sections 
for e + e~ — ► J/ipcc and e + e~ —>■ J/ipgg and find that their ratio is about 1:1.5 
at ^/s ~ 10 GeV. This result is in conflict with the current measurements by 
the Belle Collaboration, which finds that the process e + e~ — > J/ificc accounts 
for about 2/3 of all the prompt J/ifi's. We show that the discrepancy in the 
total rate as well as in the J/if> momentum distributions can be resolved by 
considering a large renormalization K factor (K ~ 4) for the J/ ipcc cross section 
and by taking into account collinear suppression in the end-point energy region 
of J/if>gg production. Detailed studies of the model predictions in terms of the 
density matrix are performed and various momentum and angular distributions 
are presented as functions of the K factors. These distributions can be used to 
determine the normalizations of each subprocess provided that the production 
and decay angular distributions do not alter much by higher order corrections. 
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1 Introduction 



As one of the simplest processes to investigate both perturbative and nonperturbative prop- 
erties of quantum chromodynamics (QCD), charmonium production at various collision 
processes has stimulated a lot of interesting theoretical and experimental works. One of 
such developments in the past few years, called nonrelativistic quantum chromodynamics 
(NRQCD) £Q which generalizes and improves the conventional color-singlet model (CSM), 
has provided a successful explanation of the CDF measurements of prompt J/ip and ip' pro- 
duction at the Tevatron |2j. Within the framework of NRQCD, the puzzle of ip productions 
in excess of the CSM prediction can be solved by introducing significant contributions from 
color-octet terms, which correspond to a gluon forming a cc pair in a color-octet state at short 
distances and then evolving at long distances into a color-singlet state along with other light 
quarks. To confirm the validity of the color-octet mechanism (COM), Braaten and Chen 
first suggested that the inclusive ip productions through e + e~ annihilation may provide an 
opportunity to observe the color-octet contributions to the cross section and angular distri- 
butions 3 j. Subsequently many detailed calculations have been performed in the literature 

H UU 13 El 13 HD|. 

Recently, BaBar and Belle Collaborations have published their experimental data for 
prompt J j ip productions fT| El • Both measurements for the inclusive processes are 
dramatically larger than the leading-order prediction of the CSM. However, there is no 
obvious evidence in support of the existence of the color-octet state reported by Belle mea- 
surement, especially in the upper end-point region of J/ip momentum distributions where 
the significant color-octet signal is predicted by previous analysis. Although higher order 
corrections may significantly soften the hard J/ijj momentum spectrum of the COM ^3], it 
is the following Belle observation 

a(e + e" -> J/tpcc) /a(e + e~ -> J/ipX) = 0.59j£ijj ± 0.12 = 0.59 ± 0.18, (1) 

which prompts us to reconsider the normalization of the J/i/jcc production in CSM. Even 
higher ratio, 

a(e + e~ -> J/^cc)/a(e + e" -> J/if)X) = 0.67 ± 0.12, (2) 

has been reported by Belle as a preliminary result based on 86.7fb _1 data set [T^]. Since it 
is unlikely that COM gives the ratio larger than a quarter, we may conclude that the COM 
contribution is sub-dominant in the process e + e~ — > J/ipX at y/s = 10 GeV. 

In this report, we assume that the COM contribution is negligible for prompt J/ip pro- 
ductions, and study if the CSM predictions for e + e~ — > J/ipcc and e + e _ — > J/ipgg can be 
made consistent by introducing the renormalization K factors for the total cross sections, 
and by considering the softening of the J / ip momentum spectrum in the J/ ipgg process. We 
give predictions for the J/ ip production and decay angular distributions, which will be useful 
to test our assumptions and to measure the K factors and the momentum distributions of 
individual subprocesses in future precision experiments. 
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We note here some other evidences that may support our assumption indirectly. The 
exclusive cross sections for the double-charmonium productions such as e + e~ — > J/iprj c , mea- 
sured by Belle, is also one order of magnitude larger than the leading-order prediction |13j . 
Since in the exclusive processes, the color-octet contribution is negligible, a large K factor 
for the e + e~ — ► J/tprj c production amplitude in the CSM may be necessary to explain the 
experimental data, as has been illustrated in Ref. [T5] . 

Lately some authors have calculated electromagnetic contributions and the relativistic 
correction [T7j, and also considered the possibility of mis-detecting the QED process e + e~ — > 
J/ijjJ/ijj at Belle |18[ITH]. However, since the considered subprocesses have the same order of 
magnitude of the original tree-level cross section in CSM, they are not sufficient to provide 
an explanation of the large enhancement of the J/i/jcc cross section. In the recent report [20J, 
Belle collaboration showed that they are capable of distinguishing J / ip from rj c and that they 
observe no evidence of e + e~ — > J/ip J/ip yet. 

Another evidence of possible suppression of COM comes from radiative T decays. At 
the amplitude level, T — > 7X is quite similar to inclusive J/ip productions except for the 
heavy quark masses. Recently, radiative T decays observe by CLEO [21] has been studied in 
the framework of the combination of NRQCD and the soft-collinear effective theory (SCET) 
|2*2*1 12*3*] . The result shows that a good fit in the end-point region of the photon energy 
spectrum is obtained only when the color-octet matrix elements are set to zero. The same 
may or may not apply for prompt J/ip production in e + e~ collisions, but we can take into 
account the collinear suppression effect at the end-point region in e + e~ — > J/ipgg process, 
which is very similar to that of the T — > ^gg in the CSM. 

In Section 2, we give the leading-order calculations of the processes e + e~ — > J/ipcc and 
e + e~ — > J I ipgg and correct some mistakes in literatures. The production of J / ipcc through 
two virtual photons is also considered and the fractions for transversely and longitudinally 
polarized J/ip are presented. In Section 3, we extract the collinear suppression effect in 
e + e~ — > J /ipgg production from the photon momentum spectrum of T radiative decays, 
and find that the observed J / ip momentum spectrum and the J / i/jcc fraction Eq. (0) can be 
reproduced by the CSM if we introduce a large K factor of K ~ 4 for the J/ipcc cross sec- 
tion while the J / ipgg process does not need a K factor significantly different from unity. In 
Section 4, we present various distributions for the inclusive J/ip productions in the different 
momentum regions and study the sensitivity of the K factor in terms of several production 
and decay angular asymmetries. In Section 5, We discuss the connection between the in- 
clusive process and the exclusive double-charmonium production. Finally, our conclusions 
are given in Section 6. In Appendix A, we present technical details of the scheme which we 
developed to calculate the helicity amplitudes of NRQCD processes by using the HELAS 
codes 

2 The Leading-order Calculations 

In the color-singlet picture, three production modes are involved in the inclusive J/ip pro- 
duction processes, as shown in Fig. 1 and Fig. 2. We refer them to the QCD J/ijjcc produc- 
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tion (Fig. 1(a)), the QED J/ipcc production (Fig. l(b),(c),(d)) and the J/ipgg production 
(Fig. 2). Other QED processes such as e + e~ — > J/if)j* where the 7* forms a lepton or quark 
pair (except for the charm pair) have been regarded as the background and removed exper- 
imentally The J/ipgg production was initially investigated twenty years ago and was 
taken as the most important production mode in comparison with the color-octet production 
in the old-fashioned "color-evaporation model" [26J. The calculations on the QCD J/ipcc 
production have been carried out in Refs. [U 03 E] and shows that the cross section is com- 
paratively smaller than that of J/ipgg production. In the case of the t-channel QED J/ipcc 
production (Fig. 1(b)), Ref. [Oj shows that the cross section is only 1/6 of the QCD one and 
the modest effects are expected. Here in spite of a suppression factor a 2 EM /a 2 compared to 
the QCD J/ipcc production, the t-channel diagrams of the QED J/ipcc process should be 
taken into account due to the enhancement of powers of y / s/2m c , which is studied in detail 
in Refs. According to the Furry's theorem, the density matrix of the interferences 

between s- and t- channel J/ipcc processes (including the s- channel QCD diagrams) vanishes 
when cc angular distributions in the cc rest frame are integrated out. We also calculate the 
corrections from s-channel QED productions (the detailed procedures of our calculations will 
be shown later) 

-3.0% + 3.5 x 1(T 3 , 

2.0% + 1.0 x 10~ 4 , (3) 

where the symbols |a| 2 , |c| 2 , \d\ 2 represent the contributions to the total cross section from 
diagrams Fig. 1(a), (c) and (d), respectively. The first values on the right-hand side of Eq. (jH)) 
represent the corrections from the interference terms 2Re(ac*)/|a| 2 and 2Re(ad*)/|a| 2 , while 
the second values represent the direct corrections |c| 2 /|a| 2 and |ci| 2 /|a| 2 . We note that the 
interference between the diagrams (a) and (c) is negative, which is consistent with the QED 
corrections to the exclusive process e + e~ — > J/ipr/ c [T7]. As expected, the sum of the above 
corrections from the s-channel QED processes contributes to the total J/ipcc cross section 
destructively by about 1%, which is sufficiently small and can be neglected safely. Therefore, 
in the rest of the paper, we refer only to the t-channel QED Jipcc production of Figl. (b) 
as the QED process. 

It has been noted that the color-octet process e + e~ — > [cc]sg —> J l^g exhibits a remark- 
ably different property with all of the above three color-singlet productions. This process 
gives rise to a hard spectrum where the J/ip momentum is almost maximal because of its 
two-body final state This has been a crucial feature used in experiments to distinguish 
the color-octet contributions. As mentioned in the last section, no such signature has been 
detected by BaBar [TT] and Belle O EB] . Recent investigation reveals that the color-octet 
contribution to the J/ip spectrum can be broaden significantly by the large perturbative 
corrections and enhanced nonperturbative effects (T3]. In this sense, one cannot rule out 
the color-octet contribution from current experimental data yet. However, as concluded in 
Ref. 114], due to the unknown color-octet shape functions for the J/ip production, the dom- 
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Figure 1: Feynman diagrams for the J/ipcc productions from e + + e~ annihilation: (a) 
the QCD production, (b) the t-channel QED production, (c) and (d) the s-channel QED 
productions. 
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J/ip + 5 cross diagrams 



Figure 2: Feynman diagrams for the J/ijjgg production from e + + e annihilation. 
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inant J/ifjcc cross section must be understood accurately before extracting the color-octet 
contributions from the data. 

We now present our leading-order calculations. The traditional approach for 5*-wave 
productions is either calculating the helicity amplitudes by using a covariant projection 
formalism in the CSM, or evaluating the squared amplitudes directly with the help of the 
optical theorem in the general NRQCD factorization framework. To the lowest order in the 
power expansion in the relativistic velocity v of the heavy quark and anti-quark, the two 
methods are identical with each other. In this paper, we develop a much simpler numerical 
method for the amplitude calculations. First we apply the program MadGraph to 
generate the parton-level helicity amplitudes for the processes e + e~ — > cccc and e + e~ — > 
ccgg, and then combine a color-singlet cc pair to be a J/if> meson by using the HELAS 
subroutines [21]. We leave the detailed procedure for constructing the helicity amplitudes 
for NRQCD processes in the appendix. Finally the integration over the phase space of the 
squared amplitudes are performed by the Monte Carlo program BASES |27j. This approach 
avoids lengthy trace computations and allows us to obtain various production and decay 
angular distributions from the density matrices that can easily be obtained from the helicity 
amplitudes. 

In our numerical analysis we use the input parameter values as follows: the center-of- 
mass energy y/s = 10.58 GeV, m c = 1.5 GeV, a s {2m c ) = 0.258, and asM(2m c ) = 1/129.6. 
The radial wave function for J /if) at the origin -R(O) is measured through the leptonic decay 
width T(J/if) -> e+e"), 

9M 2 

1^(0) | 2 = p-T( J /if) -> e+e") = 0.447 GeV 3 , (4) 

16a EM 

where T(J/ip — > e + e~) = 5.26 x 10~ 6 GeV, and Mj/^ ~ 2m c . Here we adopt the -R(O) value 
obtained from the above leading-order formula, rather than from the next-to-leading-order 
one or from some potential models [2H|, in order to define the renormalization K factor 
of the J /if) production cross section unanimously. 

With the above parameters, the leading-order cross sections for the inclusive J /if) pro- 
cesses are given as 

a® CD = 0.0897 pb for the QCD J/ipcc production, 
a QED = 0.0156 pb for the QED J/ipcc production, 

<jgg = 0.162 pb for the J/ipgg production. (5) 

The ratio for the configurations of J/tf>cc and J/if>gg, i.e.(a^, + a ?c ED ) l a ggi ls about 1:1.5, 
while cr^ c CD /<Jgg = 1 : 1.8, which are much larger than the previous results given in Refs. (3J 
El El - For cr£ CD , choosing the same input values used in Refs. jUEHH], our result agrees with 
Ref. jH] but disagrees with that in Ref. |H] and is three times larger than that in Ref. jlj. 
This missing factor 3 was also pointed out in Ref. jB]. In the case of a gg , our result is smaller 
than those in Refs. jUE] by about a factor of 2, but is consistent with that in Refs. [2H]- in 
addition, our result for a^ c ED agrees with that in Ref. 0. Both a^ D and a gg contribute 
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to the total inclusive cross section at 0(a^a% M ), while a^ c ED contributes at 0(a% M ). The 
QED cross section is enhanced by powers of \fs/2m c . As we shall see, the contribution 
from the QED J /ipcc production affects the angular distributions significantly at larger J/if> 
momenta. 

We also examine the consistency of our numerical results with the fragmentation approx- 
imation in the high energy limit; \fs ^> 2m c . In Ref. |2H], the charm quark fragmentation 
function into J /if) has been defined as 



*°° 2 <r(e+e~ — > cc) dx 



V c ^j/^{x) = hm - _ ( ^^_ — — - — JZ—( e+e J/ipcc) (6) 



and is found to be 



„ o N2 | J R(0)| 2 x(l-x) 2 (16 - 32x + 72a; 2 -32x 3 + 5a; 4 ) 



with x = lEjj^jyfs (Ej/ij, is the J/ip energy). We confirm the above result numerically 
and find that the deviation between the fragmentation formula (J7J) and the exact differential 
cross section decreases as m c / yfs at high energies. The difference is about 6.2% at y/s = 50 
GeV and 2.4% at = 100 GeV. 

In Fig. 3, we show the momentum and angular distributions for the three types of 
J /if) productions. Throughout the paper, we adopt a dimensionless variable z for the J /if? 
momentum distributions. The relation between z and the J /if} momentum Pj/^p is given 
by z = Pj/^/Pjj^, where Pff denotes the maximum value of the J /if) momentum in the 
inclusive J /if) processes, namely, PJj^ = (s — Mj,^)/2y/s ~ 4.86 GeV. For e + e _ — > J/if)cc, 



Pj/tP cannot achieve -PJyyf due to the kinematic constraint. The largest Pj/^ is ys/4 — 4m 2 
which is about 4.36 GeV, corresponding to z ~ 0.90 at yfs = 10.58 GeV. In the case of 
angular distributions, we define 9 as the opening angle between the produced J/if> momentum 
and the electron beam. Both z and 9 are defined in the e + e~ collision center-of-mass frame. 

From Fig. 3(a) and (b), we observe that although the QED J/if>cc cross section is much 
smaller than the QCD one, its effects in the high momentum region and at large | cos#| cannot 
be neglected. This is due to the characteristics of the photon fragmentation diagrams where 
one virtual photon fragments into a charmonium and another one evolves into a cc pair (see 
the left-most Feynman diagram in Fig. 1(b)). The photon fragmentation configurations 
make up 71% of the QED cross section, while the rest comes from the interference terms 
and non-fragmentation terms. In the end-point region 0.9 < z < 1, only the J/tpgg mode is 
allowed. 

The differential cross sections for the QCD J/ipcc and the J/ipgg productions are re- 
stricted by unitarity, parity, and angular momentum considerations and can be parameter- 
ized to a simple form [3] 

S(z)[l + a(z)cos 2 9], (8) 



dzd cos 9 
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where the angular coefficient a(z) is generally limited in the interval — 1 < a(z) < 1. The 
angular distributions of the QCD J/ipcc process (denoted by solid curves) in Fig. 3(b) gives 
a positive a (a ~ 0.46), corresponding to the dominance of transverse J/if> mesons, while 
the flat shape of the J/ipgg angular distribution (dotted lines) gives a ~ 0, corresponding 
to a large fraction of the longitudinally polarized J/ip [ol/^t ~ 1.8). These features can be 
confirmed in Fig. 3(c) and (d), where we show the fraction of transversally polarized J/if> 
as function of z and cos 8, respectively. In contrast to the above cosine-square behavior in 
Eq. (|SJ), the differential cross section for the QED J/ipcc production has strong enhancement 
near | cos#| ~ 1 because of the t-channel electron exchange amplitude, Fig. 1(b). 

The fraction of the transversely polarized J /if) is displayed in Fig. 3(c) and (d) where the 
transverse cross sections have been normalized by the respective total cross sections for the 
three production modes. For the J/if)cc mode, both the QCD and QED cross sections are 
dominated by the transverse J /if), and the transverse J /if) fraction increases with the J/if> 
momentum. On the other hand, an opposite behavior is obtained for J/ipgg. We point out 
here that the transverse fraction for the QCD J/ipcc production approaches to unity as z 
reaches 0.90, the maximum for the process. This feature is not clearly seen in Fig. 3(c) due 
to the finite bin size, but has been checked numerically. It indicates that all J /if) mesons 
are transversely polarized if the associated cc has vanishing relative momentum. This is 
consistent with the exclusive double-charmonium production e + e~ — > J/ipr] c where all J /if) 
mesons are transversely polarized. The scattering-angle distributions are also quite distinct 
from each other; as shown in Fig. 3(d). In the region of the small scattering angles, the 
transverse components for the QCD J/ipcc and J/ipgg processes are comparatively larger 
than those in the large-angle region. For the QED J/ipcc process, the transverse fraction 
rapidly falls off in the region | cos^l > 0.8. 

3 The Collinear Suppression and K Factors 

We note that although the contribution from the J/ipgg process is significant in the end- 
point region 0.90 < z < 1 in our leading-order calculations, there is no such apparent signal 
observed by experiments. The probable reason of this discrepancy is that NRQCD does not 
include collinear degrees of freedom at large z and therefore both the perturbative expansion 
and the operator product expansion (OPE) break down |3U]. The problem can be cured by 
combining NRQCD for the heavy degrees of freedom with the soft-collinear effective theory 
(SCET) for the light degrees of freedom |31j . A similar problem occurs in the radiative decay 
T — ► jgg at large 7 energies as mentioned in Sec. 1. In Refs. it has been shown that 
the prediction of SCET in the end-point region of the 7 energy distribution is much closer 
to the CLEO data than the leading-order calculations in NRQCD. 

In this paper instead of performing the complicated calculations within SCET, we adopt 
a phenomenological approach to obtain an appropriate end-point spectrum for the J/ipgg 
process by noting that in SCET, the operator of 7* — > J/ipgg is in the same form as that 
of T — > ^ygg. It is further observed that, the collinear Sudakov factor of T decay arises 
solely from the gluon jet function. In the leading power expansion of SCET, the jet function 
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Figure 3: Momentum and angular distributions of e + e~ — > J/ipcc and J/ipgg processes: (a) 
and (b) are for the total cross section, and (c) and (d) are for the fraction of the transversely 
polarized J/ip production. The solid, dashed and dotted lines correspond to the QCD J/ipcc, 
QED J/ipcc, and J/tpgg production, respectively. 
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depends only on the large light-cone momentum component. It is easy to find that, the large 
momentum component of two-gluon jet of T — > 'jgg decay (Mx) is roughly equal to that 
of 7* — > J/4>gg process ( (s — Mj^)/y/s). This implies that the collinear Sudakov factor 
of these two processes are quite similar. Therefore it should be a reasonable approximation 
that these two processes have the same collinear suppression function. 

Accordingly, we disentangle the collinear suppression function F(z') from the experimen- 
tal data and the tree level prediction for T — ■> jgg as * 

where T exp and Ylo are the decay width for the experimental measurement and the leading- 
order calculation respectively, z' is the normalized photon energy and can be expressed as a 
function of the invariant mass of the two gluons 



2£v . Mj 



We assume that the function F(z') can be parameterized as 

,1 for < z' < 0.5 

F(z') = { ~ ~ (11) 

;i - z')exp( Cl z' + c 2 z' 2 ) for 0.5 < z' < 1 

with two free parameters ci and C2. We find a good fit to the data and by using the LO 
calculation presented in Ref. [SB], for 

ci = 0.96, c 2 = 0.69. (12) 

The suppression function in Eq. (|llj) shows that the collinear effects in the region < z 1 < 
0.5 are negligible, the large suppression takes place in the end-point region 0.7 < z' < 1. With 
Eqs. (JUJ and (JEJ), the suppression effects in the T radiative decay can be well described. 
Finally, we obtain the modified differential cross section for e + e~ — ► J/ipgg 



^-Jj^^y^., (I3 , 

where we replace the photon energy fraction z' by the J/ip energy fraction 2Ej/^/y/s = 
J ' z 2 + AMj^/s. We expect that this is a good approximation because the energy scale of 
the two processes is roughly the same. Integrating over z, the total cross section is reduced 
to 

a gg = 0.123 pb, (14) 



*We obtain the experimental data and the tree level prediction from Fig. 3 in the first reference in 
Refs. B3|. 
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Figure 4: Momentum (a) and angular (b) distributions for the J/ipgg cross section. The 
dashed line is for the leading order calculation, while the solid line is for the calculation with 
the collinear suppression. 

which is about 24% smaller than the original value in Eq. (j5J). 

The end-point suppression behavior is plotted in Fig. 4(a) by the solid line, while the 
result of the leading order calculation is shown by the dashed line for comparison. In terms 
of the momentum fraction variable z, the suppression effect starts at around z = 0.4. As z is 
close to unity, the differential cross section drops to zero, which diminishes the disagreement 
with the Belle data. The angular distribution is modified slightly, as shown in Fig. 4(b). 
The cross section in the central region (cos 9 ~ 0) is suppressed most because that the large z 
component dominated by the longitudinally polarized J/ijj mesons is suppressed significantly. 
After the collinear suppression factor is introduced, there is no apparent disagreement be- 
tween the Belle data at large z and the NRQCD production of the J/ipgg process, without 
introducing a further overall suppression factor. 

The NRQCD prediction for the direct J/i/jcc productions is however significantly smaller 
than the Belle data. The prompt J/ipcc production cross section obtained by Belle is 
^prompt = (0. 87_to]i9 ± 0.17) pb based on D* + and D° measurements [T3]. Because the 
prompt J/ijj data contain contributions from the ip(2S) — > J/tp transitions, we estimate the 
direct J/ipcc production cross section as follows. Working in the framework of NRQCD, 
we find that direct ip(2S)cc and J/ipcc cross sections are proportional to the meson wave 
functions at the origin, and to the leptonic decay widths of ip(2S) and J/ip. According to 
Eq. (JJJ), we have 

a dir (i>(2S)cc) : a dir (J/^cc) = M} {2S) T^(2S) -> e+e") : M 2 J/f T(J/^ -> e+e~). (15) 
By using T(ip(2S) — > e + e~) = 2.15 x 10~ 6 GeV and the branching ratio for the ip{2S) — > 
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J/ipX transition fraction B = 55.7% (S3, we obtain the direct J/ipcc cross section 

&dir{J / IpCc) = (T prom pt 



1 + B{ijj{2S) -> J/^X) x ^ (25) V V 7 



Mj^riJ/l* - e+e-) 
~ 0.66 ± 0.26 pb. (16) 

Here we estimate the error of a^ r from the statistical and systematic errors of <J prornpt . 
Comparing with the result given in Eq. (J5J), a cc (J/ipcc) = a®. CD + cr® c ED = 0.11 pb, there 
is still a large gap between the experimental central value and the theoretical prediction. 
Moreover, the rate of <J cc /(a cc + a gg ) is about 0.46 evaluated from Eqs. (jSJ) and (JH|), much 
smaller than the current experimental measurement 0.67 ± 0.12 in Eq. (J2J) as indicated in 
Sec. 1. 

According to the argument given in the first section, the possible source of the disagree- 
ment may arise from the large higher-order corrections to the QCD J /ipcc production process. 
For instance, if we set as K factor to unity for both QED J/ipcc and J/ipgg processes, the 
direct J/ipcc cross section in Eq. (|TH|) requires a K factor for the QCD J /ipcc production to 
be 

K = 7.2 ± 2.9, (17) 

while the fraction of the J j ipcc process in Eq. (j2J) gives rise to 

K = 2.6+j^. (18) 

To fit to both the cross section and the fraction, we set the K factor to a moderate value 4 
for the QCD J /ipcc process, and then find 

a cc = Kaf c CD + af c ED = 0.374 pb, a cc /(a cc + a gg ) = 75%. (19) 

A part of the large K factor may come from the uncertainty in the wave function at the 
origin. By replacing the leading-order formula of the leptonic J/ip decay in Eq. (0J) by 
the next-to-leading-order one, the value of the wave function |i?(0)| 2 is changed to be 0.792 
GeV 3 [17] . while the wave function from some potential models is even larger, such as 0.810 
GeV 3 |2B| • Because the cross section is proportional to |i?(0)| 2 , this effect gives the K factor 
of 792/447 ~ 1.77 and 810/447 ~ 1.81 respectively. Another factor of two from the higher- 
order corrections to the hard scattering part may lead to the combined effect as large as 
K ~ 4. For the QED J '/ipcc case, since the fragmentation configurations are dominant, the 
ratio cx£- ED /r( J/ip — > e + e~) is approximately a constant to any order of QCD perturbation 
theory. For this reason we choose the corresponding K factor to be equal to one. The J/ ipgg 
cross section can also be enlarged by a factor of two caused by the the next-to-leading-order 
value of |i?(0)| 2 , but the higher-order QCD corrections of the hard part is unknown. For 
simplification and consistency with experimental data we set its K factor to be about one. 

We emphasize here that only the qualitative aspects of the above discussions on the K 
factors should be taken seriously. The K factors should be calculated explicitly, and the 
predictions should be confronted against experiments. In the next section we show how to 
distinguish the contributing subprocesses experimentally. 
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4 Spin Density Matrix and Sensitivity to K Factors 



In this section, we give more detailed analyses on the elements of the spin density matrix 
and various distributions. 

Since the J /if mesons are reconstructed experimentally by using the leptonic decays 
J /if) — > , e + e~, we describe the decay processes in the J /if) rest frame by the polar and 
the azimuthal angle, 9* and <f>*. The polar axis is chosen along the J/if> momentum in the 
e + e~ cm. frame. The reduced helicity amplitudes for e + e~ — > J/ipX are denoted by M\, 
where A = +1, 0, —1 is the J /if) helicity. Other particle's helicities are unobserved (integrated 
out) here. Then the production density matrix can be expressed as R\y = SM A M A *, where 
the summation is over the other particle helicities and momenta. For e + e~ — > J/if)X — > 
l + l~X, the differential cross section can be expressed as 



da- 



totai 



dzd cos 9 d cos 9*c 



= |- |(i? ++ + R—)^(l + cos 2 6*) + R o\(l - cos 2 9*) 

+ 2Re(R +0 - R o-)— ^ sin 20* cos0* - 2Im(R +0 - R -)^= sin 29* sintp* 

+ 2Re(R + _)^sin 2 rcos20*-2Im(R + _)^sin 2 ^sin20*|. (20) 

Here R\y are the functions of z and cos 9 and we use the normalization da total / dzd cos 9 = 

(R ++ +R +i?oo) after integrating over cos 9* and (f>*. From Eq. (f2T)|) . one can find that there 

are six independent combinations of the density matrix elements for the J /if) productions 
that can be measured through the decay-lepton angular distributions. 

According to the above common form, we calculate the different density matrix elements 

for all J /if) inclusive processes and find that two of them Im(R +0 — Ro-) and Im(R^ ) vanish 

after integration over the unobserved internal momenta of cc or gg system because of CP 

invariance. Among the surviving four terms, (R ++ + R ) is actually the cross section for the 

transversely polarized J /if) and i?oo is for the longitudinally polarized J /if). After integrating 
over z and cos#, three of them survive 

at = J Roodzd cos#, 

<tt = J '(R++ + R )dzd cos 9, 

: J 2Re(R + _)dzdcos9. (21) 



where a = a L + a<r is the total cross section. We also note that is from the nondiagonal 
spin density matrix elements and should be measured through the azimuthal angle of the 
decay leptons. The term proportional to Re(R +0 — Ro-) vanished after integration over cos 9. 

We plot the momentum and angular distributions for a, ax and cr + _ respectively in Fig. 
5. The solid lines are for the sum of QCD and QED e + e~ — > J/ipcc processes, where we 
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set K = 4 factor for the QCD process. The dashed lines represent the e + e~ — > J/ipgg 
contribution where the collinear suppression factor has been included. The dotted lines 
show their sum, the total J ftp inclusive cross sections. Due to its large K factor, the 
J/ipcc contribution is now dominant. Therefore, as shown in Fig. 5(a), our prediction of 
the total differential cross section versus the J/tp momentum is roughly consistent with 
the experimental observation. The small but finite contribution from the J/ipgg process 
at z > 0.9 should eventually be observed by experiments, and its magnitude will constrain 
the K factor for the J/ipgg process for which we set K — 1 in our calculation. From Fig. 
5(b), we see that the J /ipcc mode prefers large | cos#|, while the J/ipgg mode gives rather 
flat cos 6 distribution. Comparing with the transverse fractions for the QCD J/ipcc process 
displayed in Fig. 3(c) and (d), we find those for the summed J/ipcc mode in Fig. 5(c) and 
(d) are changed only slightly because of the large K factor for the QCD process. In the 
region |cos#| > 0.9, the small suppression from the QED J/ipcc component is observed. 
The nondiagonal element <r + _ is quite different for J/ipcc and J/ipgg, in particular at large 
z where <r + _ has the opposite signs. Since <r + _ is the coefficient of cos 20* depicted in Eq. 
(120)1 . it should approach zero as the direction of the J/ip momentum is along the beamline, 
i.e. | cos^l = 1. It can be seen from the angular distributions for a H plotted in Fig. 5(f). 

In Fig. 6, we extend our above analyses on the angular distributions for different mo- 
mentum regions: < z < 0.4(l-a,b,c), 0.4 < z < 0.7(2-a,b,c) and 0.7 < z < l(3-a,b,c), 
corresponding to < P J/V , < 1.94 GeV, 1.94 GeV < P JH> < 3.40 GeV and 3.40 GeV < 
Pj/ip < 4.86 GeV, respectively. The angular distribution of the total cross section da jd cos 
is given in the figures (l-a,2-a,3-a), where the J/ipgg contribution is almost flat at all z re- 
gion, and the contribution from the J /ipcc process is more pronounced at large | cos#|. The 
ratio of da/dcos8 at | cos#| = 1 and cos# = is about 1.3 at small z(l-a), 1.5 at medium 
z(2-a), and 1.7 at large z(3-a) for the J/ipcc contribution. The fractions of the transverse 
cross section dax/d cos 9 and the total distribution are shown in the figures (l-b,2-b,3-b). 
The J/ipcc contribution has generally large ctt/ct ratio, greater than about 0.7, at all cos# 
and z regions, while the J /ipgg contribution gives smaller ratio, especially at large angles and 
large z. It should also be noted that the cos 6 dependence of the o^/cr ratio for J/-0cc process 
is sensitive to the K factor for the QCD process, because the QCD and QED subprocesses 

have very different cos# distribution, as shown in Fig. 3(d). The fractions of /dcos8 

are displayed in the figures (l-c,2-c,3-c). Similar with the distributions in Fig. 5(f), the 
fraction of cr + _/rf cos 9 is zero at | cos 6*| = 1. For the J/ipgg contribution, the maximum 
value of the fraction of cr + _/dcos6 appears at cos# = and is about 0.08 at small z(l-c), 
0.04 at medium z(2-c), and 0.03 at large z(3-c). For the J/ipcc contribution, the fraction of 
a + _/dcos8 at cos# = is about 0.22, 0.12 and -0.02 at small z(l-c), medium z(2-c) and 
large z(3-c) respectively. At large z(3-c), the fraction of a + -/dcos8 is negative due to the 
large contributions from the QED J/ipcc process. Summing up, the angular distributions 
for the J/i/jcc and J/ipgg are significantly different in the various momentum regions. 

The results of Fig. 5 and Fig. 6 are shown for K = 4 for the QCD J/ipcc production. 
Since the actual value of K is still uncertain (K = 7.2±2.9 from Gdi r (J '/il'cc), and K = 2.6 jlf ;f 
from a(J/tpcc)/a(J/ipX) as shown in Eqs. ()17)) and (J18j0 . we study the sensitivity of various 
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Figure 5: Momentum and angular distributions: (a) and (b) are for the total cross section, 
and (c) and (d) are for the the fraction of the transversely polarized J/ip, (e) and (f) are for 
the fraction of cr + _. The solid, dashed and dotted lines correspond to the J/ipcc, J/ipgg and 
total inclusive productions respectively. 
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Figure 6: Angular distributions: (l-a,b,c), (2-a,b,c) and (3-a,b,c) are plotted in the intervals 
< z < 0.4, 0.4 < z < 0.7 and 0.7 < z < 1 respectively. (1,2,3-a) are for the total cross 
section, (1,2,3-b) are for the fraction of the transversely polarized J ftp and (1,2,3-c) are for 
the fraction of cr + _. The solid, dashed and dotted lines correspond to the J/ipcc, J/ipgg and 
total inclusive productions respectively. 
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distributions on K. In order to give a more direct comparison to the experimental data, we 
introduce a few kinematic asymmetries. The first one is the observable a defined in Eq. (JHJ), 
which can be obtained from the J/tp scattering angle distributions. Although the angular 
distribution for the QED J / ipcc production is not a standard 1 + a cos 2 9 function, for 
comparing with experimental data, we can still use an effective function just like Eq. (jHJ) 
to fit the curves for the QED mode in the central region, |cos# < 0.9|, and then obtain 
an effective S®, ED and a^ ED ^ . The angular distribution for the QED mode is described 
effectively by 

^ QED 

^ + eWel- (22) 

Due to the small cross section of the QED J/ipcc and according to the traditional experi- 
mental data analysis, the technique is sufficient to give a good accuracy in comparison with 
the experimental data. 

Simply summing over Eqs. (jHJ) and (|22|) . we work out the effective angular coefficients 
for J/-0CC production and the total inclusive process respectively, 



QQED QED , oQCD QCD oQED QED , oQCD QCD , a n 

qQED i qQCD ' a ~ qQED i qQCD , q 

Occ ~T~ iJcc Occ "r iJcc ~r >Jgg 



The second observable is relative to the polar angle 9* in J /if) leptonic decays. One can 
integrate over all other parameters but leave only 9* in Eq. (|20|). then obtain an reduced 
formula 

S*[l + a*cos 2 r], (24) 



d cos 9* 



which is analogue to Eq. (jHJ). The coefficient a* has a simple relation with the total cross 
section a and the transverse cross section ctt- 

a = — . (25) 

2cr — o T 

The observables a and a* versus the K factor for the QCD J/ipcc process are displayed 
in Fig. 7 in the full momentum spectrum < z < l(l-a,l-b) as well as in the moderate and 
large z regions 0.4 < z < 0.7(2-a,2-b) and 0.7 < z < l(3-a,3-b) respectively. The integrated 
luminosity in the region < z < 0.4 at the B factories comes from the off-resonance data 
and hence is one order smaller than those in the moderate and large z regions. It is hard 
to measure the observables precisely in < z < 0.4, so we do not present them here. 
The dashed, dotted and solid lines correspond to the J/^cc, J/ipgg and the total inclusive 
productions respectively. The a versus K is given in figures (l-a,2-a,3-a). The a parameter 
for the J/ipgg contribution (independent of the K factor) is almost zero at all z region, 

^Here the fitting program PAW (Physics Analysis Workstation) and x 2 method have been utilized to 
obtain the a values. 
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which is consistent with the flat curves for the total distribution cross section in Figs. 5(b) 
and 6(1 ). For the J/ipcc contribution, a is sensitive to the K values, in particular 

at small K and at large z(3-a), it varies from 0.9 to 0.6 for 1 < K < 5. This is due to the 
large QED J/ipcc contribution in large z region. The a parameter for the total inclusive 
process varies from 0.25 to 0.4 for < z < l(l-a), from 0.3 to 0.42 for 0.4 < z < 0.7(2-a), 
and from 0.32 to 0.52 for 0.7 < z < l(3-a). a* versus K is shown in figures (l-a,2-a,3-a). 
The a* parameter for the J/ipcc contribution is about 0.4 at all z region almost independent 
of K, while for the J/ipgg contribution it is about —0.5 for 0.4 < z < 0.7 (2-b), —0.7 for 
0.7 < z < 1 (3-b). For the total inclusive J /if} process, the a* is more sensitive to the K 
factor than the a. It varies from -0.82 to 0.24 for < z < l(l-b) and 0.4 < z < 0.7(2-b), 
and from -0.7 to 0.2 for 0.7 < z < l(3-b). We also plot the Belle data with a one-sigma 
error which fall into the hatched regions. As shown in Fig. 7, our theoretical calculations do 
not satisfy the current data very well. The data for a are close to the case of J/ipcc only, 
while for a* there is a tendency towards the J/ipgg predictions. However the errors of the 
experimental data are still large, and more accurate measurements are needed. 

The third asymmetry can be obtained by integrating over all other parameters but leave 
only the azimuthal <f>* in Eq. (j^Uj) and therefore the off-diagonal density-matrix element o~ + _ 
can be measured. From Eq. (|20j). we find that <r + _ is proportional to sin 2 0*. We also notice 
from Fig. 6 that a + _ is much larger in the small |cos0| region. In order to enhance the 
signal as large as possible, we impose the angular cuts | cos0|, | cos9*\ < 1/4, obtaining 



da- 



totai 



cos 01, 1 cos 9* 



1 

< - 
- 4 



512vr 



(94a - 45(x T + 47(x + _ cos : 



(26) 



To retain the information of <f>*, we split the phase space in terms of <p* into eight parts and 
integrate over them separately. After recombination of the eight integrations and normal- 
ization by the total cross section, the asymmetry is constructed as 



A \ I cos 01, 1 cos0* 



1 

< - 
~ 4 



T 3 


' f(n/2A 
Jn-w/2 


-1/4)t d 




x (da tota i/d(f)*) 


T 3 


■ f (n/2H 
Jnw/2 


-1/4). d 




x {da t otai/d(f>*) 



94a, 



947TCT — 457T0Y 



(27) 



The asymmetry A is generally small except at very high z, because | cr H | for the J/ipcc 

production is quite large near z ~ 0.9 which is shown in Fig. 5(e). In Fig. 8, we plot the 
sensitivity of the K factor for the asymmetry A in the interval 0.85 < z < 1. A for the 
J/ipgg process (the dotted line) is almost zero while for the J/ipcc (the dashed line) can 
reach -0.11. Therefore A for the total inclusive production (the solid line) is sensitive to the 
K factor due to the enhancement of the J/ipcc component with the increasing K value. 
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Figure 7: Sensitivity of the K factor for a and a*: (l-a,b), (2-a,b) and (3-a,b) are plotted 
in the intervals < z < 1, 0.4 < z < 0.7 and 0.7 < z < 1 respectively. (1,2, 3-a) and 
(1,2, 3-b) are for various asymmetries a and a* respectively. The dashed, dotted and solid 
lines correspond to the J/^cc, J/ipgg and the total inclusive productions respectively. The 
experimental data with a one-sigma error fall into the hatched regions. 
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Figure 8: Sensitivity of the K factor for the asymmetry A in the interval 0.85 < z < 1 . 
The dashed, dotted and solid lines correspond to the J/i/;cc, J/ipgg and the total inclusive 
productions respectively. 

5 Relation between the Inclusive and Exclusive Pro- 
cesses 

Since the introduction of the large K factor for the QCD J/ipcc process could provide a 
reasonable account of the inclusive J/ip production data, especially the total rate and the 
J/ipcc fraction, it may also give us a hint of the solution to the puzzle of exclusive J/ipr) c 
production, where the observed cross section exceeds the NRQCD prediction by one order of 
magnitude. For instance, the J/if>rj c production cross section can be enhanced by a squared 
of the K factor for the J/ipcc production process. This can happen if the hard part receives 
large higher order corrections only when cc pair is parallel. Numerically, if we take K = 4 
for the inclusive production, the K factor for the exclusive one can be as large as 16, which 
is consistent with the estimate given in Ref. [Tfij . 

In the case that J/ip mesons take maximum momenta in e + e~ — > J/ipcc where cc pair 
have the same momentum, the inclusive process matches the exclusive process. Therefore 
we expect that near the high end point of the J/ip momentum, the inclusive process should 
have similar angular distributions with the exclusive one. Notice that Eq. ()20|) is generally 
valid for both inclusive and exclusive J/i/j processes. Similar to what has been done in the 
last section, the angular distribution for e + e~ — > J/ipr] c can be expressed as a more definite 



form PH] j 



da 



(1 + cos 2 #)(1 + cos : 



2 9*) - sin 2 ^ sin 2 9* cos 2(f)*, 



(28) 



d cos 9d cos 9*d(j)* 
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where the definitions of 6, 6* and <fi* are the same as in Eq. (|2Up. Comparing the above 
equation with Eq. (|20|1 . we get some general information: 

• It is well known that there is no longitudinally polarized J / ip produced in the exclusive 
J/iprjc process due to parity conservation; cr L /a — > in Eq. (J3TJ. The feature is also 
reflected in Eq. fl2~£J) . 

• After integrating over 9* and 0* in Eq. ()28|). the total (transverse) cross section a is 
proportional to (1 + cos 2 6). It means that the coefficient a defined in Eq. (JHJ) is equal 
to one in the exclusive limit. 

• The nondiagonal density-matrix element <x + _ is proportional to sin 2 8. 

• If the scattering angle 9 is fixed at tt/2, the rate for the diagonal and nondiagonal 
elements, <7 + _/cr, can achieve the smallest value -1. 

We confirm that the above behaviors are reproduced in the high J / ip momentum limit of the 
QCD J/ipcc production cross section. Detailed studies of the inclusive J/ipcc cross section 
at the highest J/ip momentum region may reveal the existence of an additional large QCD 
enhancement factor in the exclusive limit of the process. 

6 Conclusions 

In this paper, we have calculated the inclusive J/ip productions in e + e~ annihilation at 
the center-of-mass energy 10.58 GeV. Within the framework of NRQCD, we performed the 
leading-order calculations for the three production modes: the QCD and QED J/ipcc pro- 
cesses and the J/ipgg process. The cross section ratio for the J/ipcc and J/ipgg processes 
is found to be 1:1.5, larger than some previous estimates jU E]. In order to explain the 
Belle experimental data [T2J EH I2D] , we considered the collinear Sudakov suppression for 
J/ipgg and a large renormalization K factor for the QCD J/ipcc mode. We found that our 
predictions with the K factor of around 4 can reproduce both the inclusive J/ip momentum 
distribution and the large fraction of the J / ipcc events |I2J EH ED] • We have presented a 
complete density-matrix analysis and considered the diagonal elements as well as nondiago- 
nal ones which give rise to polar and azimuthal angle correlations in the J / ip leptonic decays. 
The various momentum and angular distributions have been plotted, and their sensitivity 
to the magnitude of the K factor has been studied. Finally we briefly discussed the relation 
between the inclusive and exclusive processes. 

Throughout this work, we have not presented any discussion on the relativistic corrections 
to the inclusive J/tp production. We know that the relativistic corrections may play a very 
important role since the charm quark is not so heavy and its velocity v is not so small within 
the charmonium. In fact, with such corrections, the cross section of the exclusive process is 
doubled 03- In the case of inclusive processes where only one charmonium is produced, the 
relativistic corrections should not exceed that for the exclusive one. Therefore if one knows 
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how large the relativistic corrections affect the total cross section, a significant fraction of 
the large K factor may be explained. 

Before closing the discussion, we note here that we find no evidence of the color-octet 
contribution to the prompt J/ip production process. If its contribution is significant, the 
color-singlet contribution to the J/ipgg process should be further suppressed in order to 
keep the large fraction of the J/ipcc process. There is a possibility that the color-octet 
processes contribute mainly to the ip' production and that its contribution to direct J/tp 
production is suppressed. Detailed studies of prompt J/ijj and ip' production data in e+e- 
B factories may reveal the nature of charmonium production dynamics. 
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Appendix A: Helicity amplitude for NRQCD 

In this appendix, we present explicitly how we calculate the helicity amplitudes for 
NRQCD processes by using the HELAS subroutines 

Consider the inclusive production of a quarkonium H with momentum P and helicity A 
in colliders, a + b — > H\(P) + X, where a and b stand for the initial particles which can be 
e + e~ or a parton pair and X denotes additional final states. If if is a S'-wave bound state, 
to the lowest order in the velocity v, both quark and anti-quark within the quarkonium are 
on-shell and carry the momentum P/2. Hence the total amplitude can be separated into 
two independent parts as follows: 

M(a + b^H x (P)+X) = ]T A(a + b^Q\ i (P/2) + Q{ 2 (P/2) + X) 

Xi,\2,i,j 

x^Q^P/2) + (ft, (P/2) - P A (P)), (A.l) 

where A is the helicity amplitude in the parton level and T is the amplitude for forming a 
bound state. Here Q(Q) should be charm or bottom quark(anti-quark) with helicity A1XA2) 
and color i(J). The momenta and helicities for the other particles have been suppressed in 
Eq. (jA.lj) . More explicitly, A and T have forms 

A = <(P/2)r<(P/2), 
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F = C(^/2)^(75«+^(P)/3)<(P/2). (A.2) 

Here T depends on various production processes and has a complicated structure which 
consists of the wave functions of initial particles a and b as well as the final particles in 
X. On the other hand T has a simple structure where a = 1,(3 = for pseudo-scalars and 
a = 0,(3 = 1 for vectors. AT is a normalization constant related to the meson wave function 
at the origin ^(0) = R(0)/\/Att. The traditional definition of HQQ vertex is given by[ 

*(0) Sjj 
2 v / 2m 



f{P + 2m). (A.3) 



where m is the heavy quark mass. The above vertex can be obtained identically by contract- 
ing u\ x {P/2) and V\ 2 (P/2) in A with T and summing over the helicities Ai and A 2 . After 
simplification of the Dirac matrices, one can find a simple relation between the normalization 
constant and the wave function 



Both of the amplitudes A and T in Eq. (|A.2|) can be expressed in a transparent manner 
by using HELAS subroutines. In particular, the program MadGraph can automatically 
generate a Fortran code for A such that we only need to write a code for T and combine 
A and T correctly. After summing over the helicities for the intermediate heavy quark and 
anti-quark as depicted in Eq. (jA.ljl . we obtain the full amplitude M\ for different polarized 
quarkoniums. This M\ is what we need for calculating the spin density matrix in Sec. 4. 

The technique can be utilized also for the quarkonium decay, for instance J/ip — > e + e~. 
The amplitude has the form 

M(J/MP)^e + e-) = J2 HJ/MP) - cjJP/2) c^(P/2)), 

x^ Ai (P/2)4 2 (P/2)^ e + e -) (A.5) 

with 

A = v{ 2 (P/2)T'ui(P/2). (A.6) 

With this amplitude, one can easily obtain the leptonic decay width. By comparing with 
the analytic formula for the leptonic decay, we verify that our normalization factor A^ in Eq. 
flA.4|) is correct. 

We can also apply the method to calculate the exclusive processes where the additional 
heavy quark pair should form another quarkonium. The amplitude is expressed as 

M{a + b^ H${P X ) + H*{P 2 )) 
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£ £ ^(0 + 6-0^/2) + gJ a (p 1 /2) + gt 3 (A/2) + gi 4 (p 2 /2)) 

Al,A2,l,j A3,A4,fe,i 

x^i(g Al (Pi/2) + gyp/2) -> ^(po) 

x^ 2 (gt 3 (P 2 /2) + g A4 (P 2 /2) - P 2 A '(P)), (A.7) 

where the production amplitude is now expressed as 

A = 4 1 (P 1 /2)r 1 <(P 1 /2)<(P 2 /2)r 2 <(P 2 /2) (A.8) 

and T\ and Ti have similar structures as those in Eq. (|A.2|) except for different momentum, 
helicity and color indices. Summation over the helicities and color indices of two pairs of c 
and c wave functions gives the helicity amplitudes for the exclusive process e + e _ — > J/i/jt) c . 
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